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We study translational correlations of the vortex center of mass positions of the 
Abrikosov flux-line lattice in superconducting samples of finite thickness L (along the 
direction of flux-lines). The Larkin correlation lengths for the center of mass mode of 
the flux lines in the presence of point and correlated disorder are computed, and we find 
that in the case of point disorder the average (i.e. center of mass) position of flux lines 
maintains positional order on length scales which scale like ~ \/L in 2-1-1 dimensions. On 
still longer length scales, however, we find using a replica Gaussian variational approach 
that center of mass correlations cross over to a power law growth of the form r^/L, 
which should be observable in superconducting thin films. 
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1. Introduction 

There has been a great interest in the properties of the Abrikosov flux-hne lattice"'^ 
(FLL) in liigh temperature superconductors^ (HTSCs) during the past eighteen 
years. This interest was motivated by the remarkable physical properties of these 
materials, such as their high superconducting critical temperatures and their 
anisotropic, layered structure, with both properties dramatically increasing the im- 
portance of thermal fluctuations. Thus, it has quickly been recognized that HTSCs 
are an excellent system for studying the combined effects of thermal fluctuations 
and quenched disorder on flux-line assemblies. These combined effects lead to a re- 
markably rich phase diagram for the flux-line system, the phenomenology of which, 
even after many years of experimental and theoretical investigations, continues to 
pose many exciting and challenging questions. ^'^ 

In contrast to the usual phase diagram of conventional, low temperature su- 
perconductors, which shows an Abrikosov flux-line lattice on the whole region 
Hci{T) < H < Hc2{T) delimited by the lower and upper critical flelds, Hci{T) 
and Hc2{T) respectively (Fig. 1), because thermal fluctuations are now stronger, 
the vortex lattice melts into a flux-line liquid over a significant region of the mag- 
netic field (ff)-temperature (T) phase diagram (Fig. 2, upper panel). 

In the presence of disorder, most theoretical studies of fiux pinning in HTSCs 
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Fig. 1. Phase diagram of the flux line system in conventional (low-temperature) superconductors. 
Below Hci{T), the superconductor is in the so-called Meissner state, characterized by perfect 
diamagnetism, while above Hc2{T) the system becomes normal (no superconductivity). 

have concentrated on trying to understand nature of translational and orientational 
order of flux-lines as a function of temperature and applied external field, and the 
implications of the possible existence of glassy phases on the noncquilibrium dy- 
namics, especially in the presence of an external current. These studies also led to 
an investigation of these questions in the more general context of elastic media in 
the presence of thermal fluctuations, pinning (random and periodic) and external 
drive, stimulated by strong connections to other mathematically related condensed 
matter systems such as charge density waves (CDWs), Wigner crystals, domain 
walls, crystal surfaces, etc. For the particular case of flux-lines in presence of dis- 
order, the phase diagram turns out to be a very rich one, with the flux-line lattice 
losing true translational long range order (LRO) and exhibiting a highly nonlinear 
flux flow resistivity. For weak disorder and low applied magnetic flelds, this is a 
phase with quasi long range order (QLRO) and no topological defects which has 
been termed Bragg glass. At higher external fields, the Bragg glass leads way to a 
vortex glass with dislocations that improve the benefit from pinning energy (Fig. 2, 
lower panel). 

On the historical level, the first and decisive step for studying the pinning prob- 
lem using statistical mechanics was done in the remarkable work of Larkin'"''^^, where 
it was recognized that pinning can be treated perturbatively inside finite size do- 
mains (the Larkin domains) of coherently pinned flux lines, whose spatial extent 
has become known as the Larkin length Rc. Beyond this length, the disorder Hamil- 
tonian has a large number of low lying metastable "ground" states, which dooms 
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Fig. 2. Phase diagram of the flux Una system in high temperature superconductors. Upper panel: 
in absence of disorder. Lower panel: in presence of disorder. 



to eventual failure all direct methods based on straightforward perturbation theory. 
Therefore, new, more sophisticated methods had to be invented to treat the pinning 
problem beyond the Larkin length Re- 

The key quantity characterizing LRO in disordered elastic media is the roughness 



C(r) = ([u(r)-u(0)P), 



(1) 



where u(r) is the displacement of a vortex from its imdistordcd position r, and 
where the angular brackets denote averaging over thermal fluctuations and the over- 
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bar denotes averaging over the quenched disorder. Within a single Larkin correlated 
volume, one can neglect the dependence of the random pinning force f (u, r) experi- 
enced by flux lines on the displacements u. Simple analysis'^ (sec below) shows that 
C(r) ~ r^^ with the so-called wandering or roughness exponent ( = (4— d)/2 {d here 
is the dimension of the lattice). Beyond the Larkin length Rc, we cannot neglect the 
dependence of f (u, r) on the displacements u. and the periodicity of the flux line sys- 
tem has also to be taken into account. Following the pioneering work of Larkin (and 
analogous works by Lee and Rice'' in the context of CDWs and by Imry and Ma® 
in the context of random ferromagnets), significant progress in understanding the 
physics on scales longer than Rc was achieved thanks to the contributions of several 
authors, such as Villain and Fernandez®, D.S. Fisher^, Nattermann^", Bouchaud, 
Mezard and Yedidia and Korshunov^^. These results were then fruitfully applied 
to the flux line lattice problem where it was shown by Giamarchi and Le Doussal^^ 
(using a variant of the replica variational method for elastic manifolds of Mezard 
and Parisi^^) that at length scales longer than Rc the flux array is characterized by 
a logarithmic growth of flux line displacements 



both in two and three dimensions, a result which has been confirmed by an inde- 
pendent, functional renormalization group (FRG) calculation.-^^ 

An issue which is not sufficiently appreciated in the literature at the present 
time is the fact that the phonon field u(r) = u(x, z) includes all the internal modes 
of the flux lines, and that positional order is in fact best characterized by the 
displacement of the center of mass (CM) positions of vortices, which represent 
the average positions of the flux-lines. Indeed, flux-line lattices are essentially two- 
dimensional, as exhibited by the fact that both direct and reciprocal lattice vectors 
of the FLL lie in two-dimensional space. This is an essential difference between 
vortex lattices and ordinary three-dimensional solids, which implies in particular 
that in order to characterize the displacement of the average position of a given 
flux-line from its equilibrium position Rj, it is enough to consider the displacement 
uoi = roi ~ of the vortex center of mass ro^. Hence, to flnd the displacement 
of the average (i.e. center of mass) position of vortices, we need to consider the 
correlator 



Although the above distinction between the general correlator C(r) of Eq. (1) and 
the CM correlator Co(r) is of no great importance and only of academic interest 
for pure (unpinned) flux line lattices in three dimensions, below we will show that 
proper handling of the CM mode of flux lines alters the asymptotic behavior of the 
translational correlation function (1) in the presence of disorder. In order to do so, 
we shall extend the replica Gaussian variational method of reference to the case 
of samples of finite thickness along the direction of fiux lines, being careful to 
separate the CM mode from the internal modes of fiux lines. 



C(r) ~ Hr/Rc), 



(2) 



Co(r) = ([uo(r)-Uo(0)]2). 



(3) 
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Our main result is that for the fluctuations of the center of mass positions of 
flux lines the logarithmic roughness of Eq. (2) should be replaced, in flnite size 

systems (and especially in thin films) by the following algebraic growth of flux-line 
fluctuations (here the vortices are directed along the z axis, and r_L = xic + yy) 

Co{r±) ~ ^ . (4) 

Indeed, the main thrust of the argument made in rcf. 13 relied on the fact that the 
elastic propagator G{q) in presence of disorder had a wavevector dependence of the 
form 

G{q) cx 1 (5) 

in d dimensions for 2 < < 4, hence the result (2). It is important to realize that 
the elastic propagator in Eq. (5) is valid, strictly speaking, for three-dimensional 
solids (although it has been derived in the context of vortex lines). If we recall that 
q"^ = + q^ (with and q^ the wavevcctors along the directions perpendicular 
and parallel to the vortices, respectively), then it is easy to see from the above result 
that the elastic propagator Go(g) for the CM (g^ = 0) mode is given by: 

Go{q±) oc 4" . (6) 

which directly leads to Eq. (4) in d = 3 dimensions. While we argue that the 
above result, Eq. (4), is in principle valid for samples of arbitrary thickness, it is 
mostly relevant to thin superconducting films with a small enough thickness 
compared to their size L± in transverse directions (as long as ^, where ^ is 

the coherence length along the direction of fiux lines, so that flux lines have both 
CM and internal fluctuations). In such films, the above equation can be of direct 
experimental relevance and will lead to the complete destruction of translational 
long range order on length scales r± > L^. 

The rest of this paper is organized as follows. In Sec. 2, we shall start our 
investigation by carefully defining the CM and internal degrees of freedom of the 
flux lines and the associated phonon fields. We shall then briefly discuss the case 
of disorder-free flux line lattices in 2 -|- 1 and 1 + 1 dimensions. In Sec. 3, we shall 
consider the effect of an external pinning potential, both perturbatively, and using 
a replica Gaussian variational approach. In Sec. 4, we briefly comment on positional 
order for the CM mode for moving flux line arrays in disorder, and in Sec. 5 we will 
present our conclusions. 

2. Elasticity in the center of mass representation 

In order to fix ideas, let us consider a d-dimensional superconducting sample in an 
applied magnetic field H = Hz. The sample thickness in the z direction will be 
denoted by L. Vortex trajectories will be parametrized by the d-dimensional vector 
ri{z) = (xj(z), z), where Xj(2;) denotes the transverse position of the ith flux line at 
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height z, and the number of transverse dimensions will be denoted by d_L = (d — 1). 
We shall introduce the following decomposition of Xi{z) in Fourier modes^^'^^ 

oo 

Ki{z)= J2 MQn)e"'-', (7) 

n=—oo 

where g„ = 2mr/L and where the Fourier components x(g„) are related to Xi{z) by 

^{qn)^j£dzMz)e"'-', (8) 
as can be verified by using the orthogonality relation: 

rL 

dze'i"''{e'i"'')* = LSn^m- (9) 

/o 

The above Fourier decomposition is similar to the decomposition of internal modes 
into Rouse modes commonly used in polymer physics^^. For the developments that 
will follow, it will prove useful to write Xi{z) in the form 

Xi{z) = Xoi + Uii{z) , (10) 

where 

xoi = Y / dz Xiiz) (11) 
^ Jo 

is the position of the center of mass of the ith flux line, while (here c.c. denotes 
complex conjugation) 

oo 

uiiW = EW9n)e'''"" + c.c.} (12) 

is the displacement of the flux-line at height z with respect to the CM position xoi. 
Note that this last quantity itself is a dynamical variable, since the CM position 
Xoi of the i-th flux line itself fluctuates around its ideal lattice position Xj. We thus 
see that the displacement of the flux line at height z with respect to its equilibrium 
position Xj is given by 

Ui{z) = Xi{z) - Xj , 

= uoi + uii(^;), (13) 

where the 2;-independent quantity Uoi = xoi — Xj is the displacement of the center 
of mass position with respect to the equilibrium position X^ (see figure 3). 

In order to be able to define mathematically tractable models, one has to perform 
a continuum limit, whereby the individual displacement fields {ui(^;)} are replaced 
by a smooth interpolating function u(x, z). The displacement field u(x, ^;) will be 
decomposed in the following Fourier modes:^° 



u(x, z) 



oo p 

^ / u(q,g„)e^('»--+«"^), (14) 
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Fig. 3. Definition of tfic center of mass and internal modes. The solid and dashed vertical lines 
indicate the location of the equilibrium and center of mass positions, respectively. U()(x) is the 
displacement of the center of mass position with respect to the ideal position x. ui (x, z) is the 
displacement of the flux line with respect to the center of mass position at height z, and u(x, z) 
is the diplacement of the flux line at height z with respect to its equilibrium position. One has 
u(x, z) = uo (x) + ui (x, z) . 



where stands for J^^ d'^^(\/ {2itY^ (note that q here is a rf_L-diinensional vector), 
and the integration is over the first Brillouin zone of the flux hne lattice. The Fourier 
components {u(q, Q'„)} can be related to u(x, 2:) through: 

u(q, qn) ^^£dzjdx u(x, z) e-'('»"^+«"^) , (15) 
as can be verified by using equation (9) and the fact that 
J dx e'('i+'»'>" = {2TTf^ (5(q + q') . 



I 



I 
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In the above representation of the elastic displacements u(x, z), the usual expres- 
sion of the elastic matrix^^'^ (here Paffio) = Qa 10 /<f and -P^(q) = <^a/3 — Pafji^) 
are the usual projection operators) 

^a/3(q, qn) = ^L(q, qn)Pa0{fl) + ^T(q, 9n)f J/3(q) , (16) 

remains unchanged, except for the appearance of discrete g„ modes (instead of a 
continuous Qz variable) and for a slight change in the expression of the longitudinal 
and transverse components which, due to our specific Fourier transform convention, 
are now given by : 

^L(q, Qn) = L{cnq'^ + cml) , (17) 
^T(q, qn) = L{cm^ + caql) . (18) 

In the above definitions, cn, cee and C44 are the compression, shear and tilt mod- 
uli of the FLL, respectively. In principle, the compression and tilt moduli have 
spatial dispersion, i.e. they are wavevector-dependent-^^. Here we shall neglect this 
dispersion altogether, as it does not affect the long distance behavior of correlation 
functions. 

The elastic Hamiltonian can now be written in the form 

Hel = 7; X] /^a(-q,-9n)^a/3(q,9n)W/3(q,an)- (19) 

It is to be noted that the elastic matrix $(q, g„ = 0) of the center of mass phonons 
in three dimensions 

$„/3(q, 0) = icii g2 pi^(q) + Lcee q^ P^^piq) (20) 

has the same form as the elastic matrix of a two-dimensional lattice of ordinary 
particles with (two-dimensional) compression and tilt moduli given by Cn = Lcn 
and Cflg = Lcqq respectively. From the expression (19) of the elastic Hamiltonian, 
we can easily write for the thermal average {ua{(l, qn)uf3{q', qm)) the following result 
(we henceforth use units such that Boltzmann's constant ks = 1) 

{ua{ci, qn)Uf){c{, qm)) = {2i^)'''^ Sn-m^^iq. + q') rGa/3(q, qn) , (21) 

with the elastic propagator 

G«/3(q,<Z„) = [^L]-'PtM) + m-^Plpiq). (22) 

If we call Go, a/3 (q) the elastic propagator of the center of mass mode Ga^(q, qn = 0), 
then it can easily be verified that the mean square relative displacement G(x) = 
([u(x) — u(0)]^) can be written as the sum 

G(x) = Go(x) + Gi(x), (23) 

where 

Go(x) = 2T^ / Go,„a(q){l - cos[q • x]} (24) 
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is the mean square displacement of the center of mass mode, while 

Ci(x,z)=2T^^ / G«„(q,g„){l -cos[q-x + g„2]} (25) 

is the corresponding quantity for the internal modes. In what follows, we shall be 

mainly interested in the center of mass relative displacement Co (x) whose behavior 
has been largely overlooked in the past, and which will turn out to control the large 
scale asymptotics of C{x,z) and the long range translational order in the vortex 
system. 

Using the properties of projection operators, one can easily verify that the trace 
of the elastic propagator Go (q) is given by : 

VGo,«a(q) = 7(^ + ^), (26) 

and hence we obtain in 2 + 1 dimensions : 



Co(x) 



- [ ^ [l-Jo(g|x|)], (27) 
>6 Jq q 



where we took into account the fact that coe «C cii (and hence that c^^ is negligibly 
small compared to c^q) for most HTSC. The large distance behavior of the integral 
on the rhs of the last equation can be obtained in a standard way^^, with the result : 



L 



^ ^ [1 - Jo(g|x|)] = In A|x| + « + 0((A|x|)-3/2) 



/o Q 

where k = J^du{l — Jq{u))/u'^ ~ —0.116 is a numerical constant. Hence, the 
correlation function Co(x) takes the form^"' 

Co(x) ~ [lnA|x| + k] . (28) 

Using the fact that cee = £o/4a^, where eq = (<?^o/47rA)^ (here 0o = hc/2e is the 
flux quantum^^ and A is the London penetration depth) and a is the vortex lattice 
spacing, we obtain 

lCo(x,^):^ (^) [lnA|x|+K], (29) 
with the characteristic temperature 

To ~ ^ Leo ■ (30) 

The fact that To is proportional to L is an indication of the fact that the CM mode 

with compression and shear moduli Cn = Lcn and Cge = Lc^e respectively, is much 
stiffer than the internal modes. The characteristic temperature Tq of equation (30) 
is in general a very large temperature for macroscopic sample thicknesses. Using 
the result^ 

(K/A) = ^-Q*^^ ^ 1.964 X lOS ,^ _ 

' ' [A(T)(A)]2 [A(0)(l)]2 ^ tJ' ^ ' 
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with A(0) = 14001 (typical of HTSC) and L = 1cm, we find tiiat Tq is of the 
order of lO^^K. Even in low temperature superconductors with much higher values 
of A(0) (typically A(0) is of order lO** to 10^ A in these materials), Tq is still very 
high for bulk samples. We therefore conclude that, even though it may appear at 
first sight from equation (29) that the center of mass positions of the vortices have 
only quasi long range order, because the prcfactor (T/Tq) is extremely small for all 
temperatures of interest on one hand, and of the very slow variation of the logarithm 
function on the other, it can be said that the center of mass mode possesses true 
long range order for all realistic experimental situations. 

The situation is, however, much less clear in 1 + 1 dimensions (i.e. flux lines 
confined to a plane). There we find 



Co{x) 



T dq 1- e*«^ 
L J_j^ 2tt ciig2 

T\x\ /-^l"! , 1-cosu 
du t; 



ttLc 
Tlx 



-J 

11 Jo 

(32) 



•^11 = — : ~ Klyi ' (^^) 



2Lcii 

where in going from the first to the second line we used the change of variables 
q = A|x|, and where we obtained the last line by sending A|x| to infinity and using 
the rcsult^^ du{l — cofiu)/v? = -k 12. Here again, the survival of positional order 
in a given sample of thickness L will depend on the numerical value of {T/ttLch). 
If the vortices are so far apart from each other that the condition a > A is satisfied, 
then one can use the following expression for the compression modulus^^'^ 

- ^T^n^ l_ 

C44 (1 — Dn) 

where n = 1/a is the density of flux lines, and £> « A is the range of the interaction 
between vortices. Assuming that the tilt modulus C44 is of order neo, and neglecting 
the factor (1 — Dn)~'^ which is of order unity for o > A, we obtain* 

cii ~ ^ . (34) 

Using this expression in equation (32), we are led to the following numerical estimate 
for the displacement correlation function: 

2^ \WV ^' ^ ^ 

where the expression (31) of £0 has been used. Thus we sec that, even at very low 
temperatures (of order, say, a few kelvins) and large sample thicknesses [L ^ 1cm = 
10^ A), long range order will be destroyed on a relatively short distance, of order 
a few lattice constants a. This is an important result in view of the fact that it 
is a generally accepted fact that line lattices in 1 + 1 dimensions have quasi-long 
range translational order*^. Here we have just shown that there might be physical 
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regimes where the very notion of an ordered lattice in the 1 + 1 geometry becomes 
questionable, even in the absence of disorder. 

For higher vortex densities, no closed form for the elastic constants cn and C44 
in a film of finite thickness, satisfying the correct boundary conditions for magnetic 
fields and supercurrents at the surface of the film^^, seems to exist in the literature. 
This case obviously deserves further attention but will not be considered in any 
detail here. It is however clear that, if Cn becomes of order neo when a <C A, then 
the center of mass correlation function will behave like ~ {T/To)\x\a, with Tq given 
by equation (30), and long range order will be effectively recovered again at all 
temperatures and realistic sample sizes. 

3. Positional order in the presence of a random pinning potential 

We now turn our attention to the more interesting problem of the positional order 
of the FLL in the presence of disorder. To start with, we shall first apply Larkin's 
argument to the center of mass mode of the flux lines, before presenting a more 
careful replica analysis in Section 3.2. 

3.1. Larkin's argument applied to the center of mass mode 

Following Larkin's original work^'^^, wc assume that the action of the random im- 
purities on the FLL can be represented by a random pinning force F(r) whose 
probability distribution is assumed to be Gaussian with zero mean and correlations 
(we remind the reader that the overbar denotes an average over the disorder) 

F„(x, ^)i=>(x', z') = W <5(x - ^)5{z - z'), (36) 

or, in Fourier space, 

Fa{ci,qn)Ffsici',qm) = LW Sc^is (27r)2(5(q + q') (37) 

The Hamiltonian of the flux line system is now given by: 

H = Hei- J^^ dz j dx F(x, z) ■ u(x, z) , (38) 

where Hd is the elastic Hamiltonian of equation (19). Note that, as it stands, the 
random pinning force F(x, z) only depends on the "ideal" positions of flux lines "x" 
and not on their displacements u(x, z). As wc will sec below (in paragraph 3.2), it is 
more appropriate to represent the random impurities by a random pinning potential 
V'(x + u(x, z)), and the random force term in equation (38) can be thought of as the 
leading term in a Taylor expansion of the random potential in terms of the phonon 
field {u(x, z)}. 

At low temperatures, the statistical mechanics associated with the Gibbs mea- 
sure exp(— iJ/T) is dominated by those configurations of the displacement field 
{u(x, z)} for which H is minimal Minimization^"^ of H with respect to 
{u(q, g„)} leads to the result 

Wa(q, Qn) = Ga/3(q, g„)^>(q, Qn), (39) 
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where G'a/3(q, is the elastic ptopagator of equation (22). It then follows that the 
correlation function {Ua{r)u0{r')) averaged over the configurations of the random 
force (36) is given by: 

(u„(r)u^(r')) = ^ / / Gaj{ci,qn)Gi3s{ci',qm) X 

X F„(q,g„)F/3(q',g„)e^('i-^+«"-) e*(q-x'+9™.') (40) 
Using expression (36) of the random force correlator, we obtain the following result 

{u^{r)u0{r'))=WLj2 f Ga-y(q, <Z„)G;3^(-q, -9n)e'''-("-"'' e'«"(^-^'), (41) 

which implies that the relative displacement correlator C{x,z) = ([u(x, z) — 
u(0, 0)]^) is given by 

C7(x, z) =2WLJ2 I Ga/3(q, an)GaM-q, -g„) [1 - e^(i-^+«"^)] . (42) 

Now if we use expression (22) of the clastic propagator G'(q, (/„), we find that the 
usual expression^'^^'^ of the correlator C(x, z) in d dimensions 

^/ ^ oTT^ f ^''^ 1 - cos(kj_ • x + fc^z) 

C(x, z) ~ 2VK / 7— —7 ^ :f (43) 

is replaced by (we remind the reader that d± = d — 1) 

C(X,.).^ V / l-C0s(q.X + ,.)_ 



n— — 00 



Obviously, if we transform the summation in this last expression into an integral we 
recover the rf-dimensional result (43), with the long distance behavior C(r) ~ r*"''. 
However, here we are mainly interested in the fluctuations of the CM mode which, 
we would like to argue, are the most relevant ones for positional order in samples of 
finite thickness (more specifically, in samples whose thickness L along the direction 
of flux lines is much smaller than their transverse size L±). We shall therefore write 

C(x,z) =Co(x) + Ci(x,z), (45) 

where 

_2W f d'^-d l-cos(q-x) 

is the correlator of the center of mass positions, while 

^ , , 2W ^ f d'^-Lq 1 - cos(q • X + g„z) 
Ci(x,z) = — )^ / — 

corresponds to the internal modes of the FLL. Taking d± = 2 and transforming 
the above sum into an integral, we obtain the usual large distance result in d = 3 
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dimensions^'^^''^ (recall that here, we neglect the dispersion of the elastic constants; 
taking the dispersion of C44 into account leads to a different behavior at short 
distances): 

'-66 ''44 

On the other hand, we have for the center of mass positions : 

^ r ^ 1- Jo(g|x|) 

<^o(x) = —2- / dq 



r ..1^, (48) 



where, in going from the first to the second line we made the change of variables 
u = qx, and where wc used the inverse of the transverse size of the system, 1/ L^, 
as an infrared cut-off to insure the convergence of the integral as g ^ 0. At large 
distances, the upper bound of the integral on the rhs of equation (48) can be replaced 
by infinity, with the result (we assume that |x| <C L±): 



du'—^. (49) 



Now, on the interval [x/L±, 1], we can approximate Jo{u) = 1 — ■u^/4 + o(u^). This 
leads to 

du—-l^^-\n{L^/x). (50) 

x/L± " ^ 

On the other hand, J^°° du ^~„a^"^ is just a numerical constant, whose value we 
shall denote by a, and which is approximately given by a ~ 0.287, so that equation 
(48) finally yields : 

Co(x)^-i-L ln{L^/x)+a' , (51) 

47r_LCgg L J 

where we defined a' = 4a ~ 1.15. Apart from the unimportant numerical constant 
a' (which can be absorbed in the cut-off R), this result is very similar to what 
has been obtained a long time ago^^"^^ for thin superconducting films. Here we see 
that this result in fact holds true for the CM mode of the fiux lines in samples of 
arbitrary thickness. 

The above perturbative analysis is correct up to the length scale Xc such that 
Cq{xc) — The quantity Xc defines a length at which the center of mass mode 
of vortices has "random-walked" a distance ^ from a given initial position. Using 
equation (51), we obtain that Xc is the solution of the following equation (henceforth 
we omit the constant a from our considerations) 



I 



I 
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A good approximation to this last quantity is obtained by replacing Xc inside the 
argument of the logarithm by a;o = (47rLcgg^^/W)^/^. We then have, to logarithmic 
accuracy: 



AttLq 



1/2 



''^'^[wHlT/xo)\ ^^^^ 
The important thing to note about Xc is that it varies with the sample thickness 

L as L^/^, and hence is very large for bulk samples. Lot us for example compare 
Xc to the in-plane characteristic length for the internal modes. For the sake of 
argument, we define as the length scale at which the equal height correlator 
Ci(x,0) = W\yi\/{2n'^c%^c\{^) reaches the value We have 

2^ 
W 

and hence 

Xc {LWf'^ 



K = ^^^^e, (54) 



(55) 



TT^^Xeo\/ln{L±/xo) ' 

where, in going from the first to the second line, we used the fact that C44 = B"^ /An 
(with B = <j)o/ci^ the magnetic induction inside the superconductor) and hence 
that^ 

ttAeo 



a/C66C44 = ■ (56) 

Now, if we use W w Sq/'^Oj where we denote by do the average distance between 
impurities in the superconducting sample, we obtain 



Xn 



Rc 7rVln(i-L/:i-o) 




(57) 



with K = A/^. Using k = 100, A(0) = 1400A, a = 500A, d = lOOA and L^ = L = 
1cm we obtain 

|| «6x 10^(1 -T/Te). (58) 

This shows that, now matter how large the internal fluctuations of the flux lines 
are, the center of mass mode of the fiux linos experiences much smaller fluctuations 
and has a much larger Larkin length than the internal modes. 

It is interesting to see how the above results are modified if correlated disorder 
is considered instead of point disorder. In that case, the random force correlations 
are given by: 



F„ (x) (x') = Wc (x - x') , (59) 

with^^ Wc ~ £0/(0^^0). An immediate consequence of the fact that the random 
force does not depend on z is that i^a(x) will couple only to the center of mass 
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mode of the flux lines, as can be easily seen using the fact that dz Ui {z) = in 
the last term of equation (38) (where we now let F(x, z) = F(x)). The Hamiltonian 
of the flux line system can now be written in the form (here and in what follows we 
concentrate on d = 3 dimensions) : 



-H"= -9")*«/3(q, 9n)w/3(q, g„) + 

+ j '«o,a(-q)*o,c</3(q)wo,/3(q) -L j dx F(x) • uo(x), (60) 

from which we obtain the interesting result (within the Larkin approximation of 
ignoring the u dependence of the random pinning force F) that correlation functions 

for the internal modes fluctuations remain unchanged and arc in fact the same as 
their pure counterparts. Only the CM mode is affected by correlated disorder in first 
order perturbation theory, and the CM part of the Hamiltonian (60) maps onto the 
Hamiltonian of a two-dimensional system with compression and bulk moduli given 
by C'li = Lcii and Cge = Lcaa respectively, subject to the random force F(x) = 
LF(x) whose variance is now given by Wc — L^Wc- The relative displacement of 
the center of mass mode Co(x) is then obtained in a very natural way as 

Co(x) = 2W, j ^ G„/3(q,0)G„/3(-q,0) [1 - e'^'^] , 
2Wc f d'^q l-cos(q-x) 



■-66 



(61) 



(27r)2 g4 

A calculation similar to the one carried out in equations (49)-(51) leads to the result 

47rci 



^o(x) = HL^/x), (62) 



which is the result obtained by Nelson and Vinokur^^. Unlike the displacement 

correlator in the presence of point disorder, equation (51), the correlator (62) does 
not depend on the sample thickness L. The resulting size Xc of the perturbative 
region for the CM mode, which verifies the implicit equation 

'^'^ = w^HlI/x,) ' ^^^^ 

does not increase with L, which translates the fact that columnar pins are much 
more effective in disrupting the long range order of the FLL than point- like pinners. 



3.2. Variational analysis 

As is well known, the Larkin analysis of the previous section breaks down beyond 
the Larkin length Xc- To find the long distance behavior of correlation functions, 

one has to resort to a more careful kind of analysis. Here we shall use the Gaus- 
sian variational method i4, 11,13,20 (QVM) to find the long distance behavior of the 
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center of mass relative displacement Co(x). The GVM has been used in the past 
to study the effect of disorder on elastic manifolds^* and flux lattices ^'i^'^o^ ^nd 
although uncontrolled (due to the absence of a small expansion parameter), it has 
been shown^'^'^^ to yield results that are in qualitative agreement with functional 
renormalization group calculations. 

Here, we shall restrict ourselves to the three dimensional case and to isotropic 
elasticity, the generalization to the more realistic case of different clastic moduli 
being relatively straightforward. Following ref.^"^, in this subsection we shall use the 
following Hamiltonian 

H = J dxdz ^c[{d^uf + (Vx • u)2] + j dxdz V{x, z)p{x, z), (64) 

where the vortex density is given by 

p{x,z) = '^52{x-Xi{z)) = ^52{x-Ri - Ui{z)) , (65) 

i i 

and where V{x, z) is a Gaussian random potential with zero mean and correlations 



V{x, z)V{x', z') = A(x -x',z- z') . (66) 
Using the replica "trick" 

'Vn 1 

lnZ= lim (67) 

n^O n 

to average the free energy over disorder, we obtain the following effective Hamilto- 
nian 

1 " /■ 

H^Sf = ^y^ / d^dz c[{d^naf + (VxUa)'] + 

- ^^'^■^ / P'^^^^ z)pb{^', ^')A(x -x',z- z') . (68) 

Now, using the following decomposition of the density^^'^^'^^ 

p{x, z) ^ po [1 - a„u„(x, z)] + po X] PK(x)e*^-" , (69) 

where pK (x:) = e~'^'"('''^) , and discarding rapidly oscillating terms, we obtain from 
equation (68) : 

= / dy^dz c[{d^viaf + {V^viaf]-^ j dxdz ^^docU'^^{x, z)d0v!'p{x, z) 

a=l a,b 

+ X ^ COS (K • (u„(x, z) - U6(x, z))) } , (70) 
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where we defined 

Ao =pI I dxA(x) , (71) 
dxA(x) e'^-^ . (72) 



The highly nonlinear form of H^ff precludes exact analysis. In order to make 
progress a common way is to use the so-called Gaussian variational method (GVM) , 
which has been developed by Mezard and Parisi^** in the context of random mani- 
folds, and was first applied to the vortex lattice by Bouchaud, Mezard and Yedidia^^ . 
This method consists in trying to find the best quadratic variational Hamiltonian Hy 
to describe the full nonlinear problem using Boguliubov's variational free energy^^'^^ 

Fy = {H-H^)y-TlnZ^, (73) 

where (• • ■)v denotes averaging with statistical weight exp{—Hy/T)/Zy, and where 
Zy = Tr(e~^"/^). As a trial Hamiltonian we take 

1 . f d'^q 

Hv = -^'^ J J^—jz{G~'^)ab{(l,qn)Ua{(l,qn)-Ub{-(i,-qn) (74) 

where {G~^)ab is an n x n matrix of variational parameters. The variational free 
energy is then given by 

= ^ E E / { W + iDSab - ^ q'] GabicL, qm) - d^TMTG)USab} 



a, 6 m ^ 
Ak 

2T 



a,b KjiO 

where we defined the difference correlation function (note that in this section, dif- 
ference correlation functions C(x, z) are defined with an extra factor l/d±, and that 
throughout the rest of this paper, no summation is implied on repeated indices) : 

Cab{^,z) = ([Ua(x,z)-Ub(0,0)]2) 

a± 

= -^E / [Gaa{ci,qn) + Gbb{(l,q„)-2cos{ci-y: + q„z)Gab{ci,qn)](JQ) 

Minimization of the variational free energy (75) with respect to [G(q, (7n)]afc for 
a^h leads to the following expression^^ 



[G ^(q,g'n)]a6 = C{q^ + ql)6ab 



Ao 



di_T 

^ exp ( - \K^C,b{^ = 0, ^ = 0)) . (77) 



q' 



Defining the self energy matrix Uab by 

[G-^ (q, qn)]ab = c{q^ + ql)Sab - Oab , (78) 
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we obtain 



An , ^ A'^Ak . 1 



Ta^b = + E ( - 2^'^«''(^ = 0,z = 0)), (79) 



while for a = b, we obtain 



(80) 

Two types of solution can exist for the variational equations (79)- (80). The first 
type is the so-called "replica-symmetric" solution which preserves the permutation 
symmetry between replicas, while the second type is a "replica-broken" solution in 
which permutation symmetry between replicas is violated. Here we shall not give 
much details on the replica symmetric solution, and refer the reader to reference^^ 
where more information can be found. The only change with respect to the lat- 
ter reference is in the form of the diagonal correlation function Gaa{<l,<ln) which, 
because the modes in our description are discrete, now takes the form 

= ^ + E^-^Ke-'--'-/^ . (81) 

In the above expression, £ is the (thermal) Lindcmann length^^ 

2^2r d'^-^d 1 

T 

, (82) 

Tica 

where the second line refers to = 2 + 1 dimensions. Due to the factor \/{q^ + (j^J^ 
on the rhs of equation (81), we find that the center of mass relative displacements 
grow as 

C0,aa(x) ~ X^-^ (83) 

which, for = 3 is nothing but the Larkin result of eq. (51). 

The analysis of quadratic perturbations, in replica space, about the replica- 
symmetric solution described above is most conveniently done using the eigen- 
value of the so-called "replicon" mode^'''^'', and such an analysis (see Appendix 
Appendix A) shows that, for the particular case d = 3, there exists a tempera- 
ture Tc = AttcL/Kq between a high-temperature, replica symmetric phase, and a 
low temperature glassy phase where replica symmetry is broken. This transition 
temperature is exactly the same as the one obtained previously in the context of 
correlated-disorder^^. Of course, for macroscopic samples {L ^ a) this transition 
temperature Tc is very large and is certainly not experimentally accessible, but 
from a conceptual point of view it is important to realize its existence for three 
dimensional flux lattices in the presence of point disorder. 

Now, in the low temperature, replica broken phase, the n ^ limit of matrices in 
replica space become functions^^'^'' of a real variable v which parametrizes pairs of 
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low lying metastable states. '^^ FoUowing^^'^^ we let G(q, g„) = Gaa{<l, Qn), z) = 
Caa{^,z), and parametrize Gab{<l,(ln) and Cab{^,z) for a 7^ 6 by G(q, g„;w) and 
C(x, z; v) respectively, with < t; < 1. The saddle point equation (79) can then be 
rewritten in the form : 

where we neglected the subdominant^^ K = component of the disorder (the 
^o<f/d±T term) and where 

C{O,0;v) = — J2 / 7^^(G(q,g„)-G(q,«„;^)). 
The algebraic rules for inversion of hierarchical matrices^^'^^ give us : 

C(0,0;.) = C(0,0;..) + 1 £ f j n ^ T^l^] U ^^^ ^''^ 

where we defined [a]{v) = va{v) — dua{u) and 



1 r rl'^~^n IT 



Taking the derivative of equation (84) with respect to v (keeping only^^ the recip- 
rocal lattice vectors K such that K = Ko),we obtain 

n=— 00 ^ ' 



a{v) f d''-^q TKl (j{v) ^-^ f d'^-iq TK^ 



L J {2^Y-^ (c<^2 + [^](„))2 L {2^Y-^ (e(^2 + g2) + ^^j(^))^ 

Transforming the sum on the rhs into an integral, we obtain 

a{v) f d'^-^q TKl ^ ^ f d'^^ TKl 

' div) 



where the ultraviolet cut-off A = 27r/a has been sent to infinity, and where we 
defined : 

_ f d'^k 1 _ (2 - Q!)7ri~''/2 

J (2^ (F + 1)2 ^ 2^^+! sin(^d/2)r(d/2) ^ ' 

In the limit L — > 00, the first term on the rhs of equation (89) will be negligibly 
small compared to the second one, and we can write 

jTKlcd^ 



1 



^iv)[^fP)Hv)r-'^/'. (91) 
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This approximation will turn out to be valid for a large range of values of the 
parameter v, covering almost the entire interval [0, 1] when L — > oo. Taking the 
derivative of equation (91) one more time with respect to v, and using the fact that 
[cr]'(t;) = va'{v), we finally obtain^^ 

where 6 = d - 2 and vq = 2TK^Cd/{{4: - d)c'^/'^). Replacing the solution (92) into 
equation (89), one can verify that the first term on the rhs of that last equation can 
indeed be neglected if 

{vo/vy^' « Lcd/c'/^Cd-i , (93) 

which gives, in d = 3 

In the limit L — > oo, we see that this last condition is satisfied nearly everywhere 
in the interval v G [0,1]. For films of finite thickness, one can easily show that 
equations (91)-(92) are still valid provided that the thickness L satisfies 

a condition which is satisfied by most thin HTSC films even for temperatures close 
to the superconducting critical temperature Tc- 

The above solution (92) for [(7](w) is a priori valid up to a critical value Vc, above 
which [a]{v) is just a constant, [a]{v) = {vc/vo)"^^^. With the knowledge of the 
analytic form of the function [cr](w), we now are in a position to find the diagonal 
part G(q, g„) such that 

2T ^ f d''-^' 



n=— OO ^ ^ 

and which is given by:^^'^^ 

c(g^ + 9„) V Jo c{q^ + Qi) + m{v)^ 
Specializing to the three-dimensional case (d = 3), if we now use the result (92) 
with = d — 2 = 1 and change the variable of integration to a; = v'^^'' = v^, we 
obtain 



"dv [a](t:) 1 



lo v^c{q^+ql) + [(T]{v) 2 V^{x + cvl{q^ + ql)) ' 
Now, using the result^^ 

dx 2 



(98) 



^/x{x + a) y/a 



Arctan(-\/ x/a) . 



I 



I 
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we finally obtain 

dv [a] (v) 



J 

Jo 



: Avct&n(vc/^Jcvi{q^ + ql)) . (99) 



/o c(g2 + ql) + [a] {v) ^ v^^T^ 

On the other hand, using the result [<t]{v) = {vc/vof' for < t; < 1 it is not difficult 
to see that 

''dv [a]{v) _ {j:-l){vjv,f ^^^^^ 



„^ c{q^ + ql) + H {v) c{q^ + ql) + {v^v^f ' 
Using the fact that \/ cVq = TKq/Attc, and approximating 

/ Attcvc \ 
Arctan , ~ 7r/2 , 

(as we are mainly interested in the (q, Qn) — * behavior of the correlation function), 
we find, in the long-wavelength limit 

which is nothing but equation (3.18) of reference^"^ in d = 3. The purpose of the cal- 
culation above, equations (84) to (101), was to show that nothing actually changes 
in Giamarchi and Le Doussal's derivation when we use the decomposition of the qz 
modes into center of mass and internal modes as wc have done throughout this pa- 
per, and that the 1/q'^ behavior of G(q, qn) in d dimension previously found by these 
authors is still valid here. From expression (101), we can easily find the expression 
of the correlation function C{x, z). We have: 

C{x,z) = Co{x)+Ci{x,z), (102) 

where 

Co(x) = —J j^-^ (1 - cos(q • X + 9„^))G(q, 0) , (103) 

is restricted to the CM mode, while 

2T . f d'^^^q 
Ci (x, 2:) = — ^ / — — (1 - cos(q • X + q^z)) G(q, g„) , 

|^(l-cos(q-x + g,2))G(q,g„) (104) 



f d^q 

~ T 

corresponds to the internal modes of the flux lines and is the correlation function 
calculated by Giamarchi and Le Doussal.^'^ We have, for the CM correlation function 
C'o(x) : 



C'o(x) = 



Att'^T f d^q l-cos(q-x) 
IjfKl J (2^0^ 



2tt r 



\q'—^. (105) 
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Changing variables from |x| to w = A|x|, we obtain 

nf^ 277 N , , /-^l^l , 1 - Jo{u) 



For large values of |x|, we can to a very good approximation extend the range of 
integration to infinity. Then, using the fact that^^ /o°° ~ Jo{u))/v?' = 1, we 
finally obtain 



C'o(x):^(-^) |x|. (106) 



27r 

On the other hand, if we use the result^^ for Ci (x) (in = 3 dimensions) 



C(x,z) = (-^) |x| + Aln|x|. (108) 



C'i(x) = -|2in|x|, (107) 

we finally obtain 

Now, the translational order correlation function at reciprocal wavevector K 

*if(x) = (p;k(x)pK(x)) (109) 

is given by^^ 

*x(x) =e-5^'^("), (110) 
which, given the result (108), implies that Cx(x) behaves for large |x| as 

TldlK^ |x| 

^0 



*K(x)~exp(^--|^ i^) |x| , (111) 



and hence we see that point disorder does destroy long range order of the FLL, 
although it does so only on asymptotic length scales |x| ~ L. The same con- 
clusion would follow from more sophisticated functional renormalization-group 
arguments^'*^'^^. 

Equations (106) and (111) are the most important results of this paper. They 
show clearly that the QLRO of the flux line lattice predicted by considering the full 
correlation function C(x, z) and integrating over all internal modes, is actually lost 
when one carefully separates out the fluctuations of the CM positions of the flux 
lines. As such, this result confirms our claim that the CM mode plays an important 
role and is in fact the most relevant one to look at when considering positional 
order of flux line lattices in samples of highly anisotropic shapes L <C L± . For such 
samples, we expect the destruction of long range range order to be observable in 
experiments. 
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4. Positional order in the moving flux line lattice in the presence 
of disorder 

The considerations of the past section can be generalized to the case of a moving FLL 

in a disordered potential. Here, wc shall restrict ourselves to the determination of 
the dynamic Larkin lengths*^ for the center of mass positions. Following standard 
arguments 46,48,42,43,44^ ^j^g ^^^^ show that the usual perturbative expression for 
the mean-square displacement of a FLL drifting with mean velocity v along the x 
direction and subject to a random force F(r) with the correlations (36); 

r'M 9W f ^''q 1 - cos(q ■ r) 
for a sample of finite thickness L should be replaced by 

In the above expressions, r] is the microscopic friction coefficient and isotropic elas- 
ticity with an elastic constant c is assumed. The mean sqiiarc displacement of the 
center of mass mode is given by the n — term in the above sum, namely 



2W r d'^-^ci 1 - cos(q ■ x) 



W = ^ 7^ /"T"„^ - (114) 



The integral in (114) has been caculated by several authors^'''*^'^^. We obtain: 

Co{x, y) = F{cx/W) , (115) 

rjLcv 

where y represents all the space variables other than x (the direction of the drive) 
and z (the direction of the magnetic field B), and where the scaling function F is 
such that F{0) = const and F{x) ~ a;^^-^''-^))/^ ~ a;(4-d)/2 ^^en \x\ > 1. Defining 
the dynamical Larkin lengths for the center of mass displacements iJg and Rq by 
the equation C{Rq,Rq) — we obtain: 

R^^rjviR^f/c. (117) 

Specializing to the case = 3, we see that Rq grows only linearly with v, in contrast 
to the transverse Larkin length of the total fluctuations^^'^^'^'^'^'* which exhibits 
exponential growth as a function of v, i?^ ~ exp(77tic^^)/A. The critical lengths 
{Rq, Rq) obtained here can however be quite large for macroscopic samples with a 
large thickness L. 



5. Conclusion 



In this article, wc have reexamined the problem of the positional order of the 
Abrikosov flux line lattice in type II superconducting samples of finite thickness 
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L, carefully separating the center of mass mode of the flux lines from their internal 
modes. While this separation turns out to be unimportant for clean systems with 

thermal fluctuations, in the presence of a random pinning potential we find that 
the Larkin length governing the growth of fluctuations of the center of mass posi- 
tions is very large for macroscopic samples and grows like vT in 2 + 1 dimensions, 
which suggests that flux line lattices on average retain translational order on length 
scales much larger than the lengths predicted so far (the latter corresponding to the 
growth of the internal modes fluctuations of the lines). Going beyond the simple 
Larkin analysis, within a Gaussian variational approximation with broken replica 
symmetry wc find that translational order of the three dimensional flux line lattice 
is destroyed, and that the logarithmic roughness predicted in previous work crosses 
over to power-law growth of the vortex displacements. Although this destruction of 
positional order only takes place on lateral length scales of order the sample thick- 
ness L, it can lead to experimentally measurable effects in superconducting thin 
films. 

Appendix A. Eigenvalue of the replicon mode and stability of the 
repUca symmetric solution 

In this appendix we consider the eigenvalue A of the so-called "replicon" mode, 
which for our problem is given by:^^'"'^^ 

A = l-;^E^'^^^"^''"^-''^^"'"^X7E/Gc(q,an), (A.1) 

where the connected propagator Gc(q, Qn) is defined by Gc(q, Qn) = Gabi^l, Qn)- 

Using the saddle point equations (79)- (80), it is not difficult to sec that G,. is given 
by Gc(q, Qn) = 1/ {cq^ + cql). Here, we shaU take Gc(q, q„) = 1/ (cg^ -|- cql + /j^), 
with a small regularizing mass fj,^, the limit )U — > being taken at the end of 
the calculation. Transforming the sum over modes (l/i)X]n ™to qz integrals in 
eqiiation (A.l) right away leads to the conclusion that in = 2 + 1 dimensions the 
eigenvalue of the replicon mode is always negative, and hence that three dimensional 
flux line lattices are always unstable to point disorder. Here, we shall instead be 
careful to separate the CM from the internal modes, upon which we obtain (here 
we only consider the case d = 2 + 1 dimensions): 

E X C.(P, . J = jj;; in (l + ^) + 5^ E (' + aTT^) • <^-'* 

m m=l ' 

and so, in the limit /U ^ we can write: 

Ge{p,Pm)l^, = ^ m (^) + G, , (A.3) 

where Gi is the constant given by 
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On the other hand, one can also easily show that 



where 



5 ^JTa^) • ^^-^^ 



Collecting all terms, we see that A can be written as 

In the limit /i ^ 0, the dominant term in the above sum is the term with K = Kq, 
which gives: 

A « 1 - _^|^e-^^o^i(A2c)-Si^-2(i-Si) . (A.8) 

Prom this last expression, we see that the sign of A when — > depends on the 
value of TK^)/4ttLc : 

• For TK^/4ttLc > 1, i.e. T > = AwLc/K^, A(/x ^ 0) = 1 > 0, the replica- 
symmetric solution is stable; 

• For TK^/AttLc < 1, or T < Tc = AttLc/K^, A(/i ^ 0) ^ -oo and the replica 
symmetric solution is unstable. 

An identical result has been obtained by Giamarchi and Le Doussal in ref.^*^ 
(where the authors adopted a separation of modes similar to ours) in the context of 
correlated disorder. Wc thus sec that the conclusions of this last reference regarding 
the stability of the replica symmetric solution result solely from the separation of 
modes into CM and internal modes, and has nothing to do with the nature of the 
(point-like or correlated) disorder considered. 
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